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Abstract. We introduce and solve a model of a thermometric measurement on a driven glassy system in 
a stationary state. We show that a thermometer with a sufficiently slow response measures a temperature 
higher than that of the environment, but that the measured temperature does not usually coincide with 
the effective temperature related to the violation of the Fluctuation-Dissipation Theorem. 
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1 Introduction 



. Thermal equilibrium is a rather subtle concept. It relies on 
' the distinction between "fast" and "slow" processes with 
, respect to a given macroscopic measurement. It follows 
' that the same system can be at equilibrium on one scale 
I and out of equilibrium on another. More strikingly it can 
I be at equilibrium but exhibiting different properties on 
' two scales at once 0]. 

The notion most intimately connected to equilibrium 

■ is temperature. It is operationally defined by the so-called 
, zeroth law of thermodynamics, which states that when 

■ two systems are in thermal equilibrium with a third one, 
then they must be in thermal equilibrium with each other. 
This allows one to define temperature as a signature of 

■ the equivalence class defined by mutual thermal equilib- 
, rium. This property makes possible the use of test systems, 

called "thermometers", to decide whether any two sys- 
tems will or will not remain in thermal equilibrium when 
brought into contact. When two systems are not in mu- 
tual equilibrium, the direction of the energy flow between 
them is determined by the second law of thermodynamics. 
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When dealing with non-equilibrium systems the chal- 
lenge is thus to produce an "effective" time-scale depen- 
dent temperature that would predict the direction of heat 
flows within this scale. 

Indeed, in the context of weak turbulence, Hohenberg 
and Shraiman H have defined an effective temperature for 
stationary non-equilibrium systems through an expression 
involving the response, the correlation and the tempera- 
ture of the heat reservoir. A closely related expression ap- 
pears in the theory of aging systems H . In a recent work 
[Q, these views have been unified for a class of out-of- 
equilibrium systems with small heat flows, which includes 
nonstationary pure relaxational systems, like glasses, and 
stationary systems, slowly driven by non-relaxational for- 
ces. The concept of effective temperatures has been further 
reviewed in refs. MM ■ In a recent experiment temperatures 
higher than the thermal bath temperatures have been ex- 
hibited in an oscillating circuit coupled to an aging glyc- 
erol sample after a quench M. 

In the present work we analyze the process of a ther- 
mometric measurement in a glassy system, by means of 
an exactly solvable model. We restrict ourselves to the 
stationary non-equilibrium regime of a driven system. We 
consider a simple system and a thermometer, both de- 
scribed by Langevin equations. By taking advantage of 
Time- Translation Invariance (TTI) the Langevin equation 
is transformed into an algebraic equation in Fourier space, 
that can be analytically solved. 
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Our thermometer is a simple physical system coupled 
to its own heat bath, which is different from the thermal 
bath of the observed system. We suppose that during the 
measuring time, the two systems are brought into contact, 
each being coupled with its own thermal bath. We then 
monitor the exchanged energy between the system and 
the thermometer in the stationary regime. The reading of 
the thermometer corresponds to the temperature of the 
heat bath of the thermometer for which the net energy 
flow between the system and the thermometer vanishes. 
We discuss the relation of the measured temperature with 
the effective temperature defined in refs. [g[Q]. 

In section 2 we recall the generalization of the Fluctua- 
tion-Dissipation Theorem to nonequilibrium systems and 
show how it defines an effective temperature. In section 
3 we describe the general measurement procedure. In sec- 
tion 4 we specify the procedure for the measurement of 
the effective temperature of an asymmetric spherical SK 
model. Finally, section 5 is devoted to the analysis of the 
results obtained for this system. 



2 Fluctuation-Dissipation Theorem and 
effective temperatures 



At equilibrium, we can use the FDT to compare this sus- 
ceptibility to the correlation function: 



X(r) 



l-g(r) 
T 



(6) 



(We are considering magnetic systems for which C(0) = 
1.) Then, a parametric plot of xi^) '"S- C(t) yields a 
straight hne with slope equal to — 1/T. 

A certain class of out of equilibrium systems with very 
slow dynamics exhibits an aging regime in which the FDT 
is violated in a very specific way. Following an initial 
quench of temperature, these systems fall out of equilib- 
rium and do not reach it again, even on macroscopic time 
scales. The longer the time i^ elapsed since the initial 
quench, the slower is the response of the system to a given 
perturbation: the system ages. This phenomenon also ap- 
pears in the time-correlation functions. In these systems, 
even in the limit iw — * cx), a parametric plot of x(i, ^w) vs. 
C{t,t^) does not yield a straight line with slope —1/T as 
in equilibrium. 

Very similar features appear in some stationary non- 
equilibrium (driven) systems [§|4||,^,H , like the one shown 
in figure |l|. Let us introduce a parameter v which measures 



According to |1, the definition of an effective temperature 
for non-equilibrium systems can be related to the viola- 
tion of the Fluctuation-Dissipation Theorem (FDT). Let 
us consider a system (described by the Hamiltonian H) 
subject to a time-dependent perturbation of the form: 



H — > H -h{t)0, 



(1) 



where O is an extensive operator. The correlation function 
C{t,t') of O is defined by 



C(i,0:-(O(t)O(i-)), 

■.^{0{t)0{t'))-{0{t)){0{t')), 



(2) 



while the corresponding response function R{t, t') is given 



Sh{t') 



(3) 



h=0 



For a system at equilibrium with a thermal reservoir at 
temperature T, Time- Translation Invariance (TTI) inti- 
mates that both the correlation and the response func- 
tions depend only on the time difference r between their 
time arguments [t = t — t') . On the other hand, the FDT 
entails a relation between the response and the correlation 
functions: 



T dt' 



T di 



Experimentally, one usually measures the (time-) integra- 
ted susceptibility: 



X(r):= / dr'Rir'). 
Jq 



(5) 



X 0.5 




Fig. 1. Plot of X vs. C of the asymmetric spherical SK model 
for Ti = 10, J — 20 and different values of the asymmetry 
parameter v. The susceptibility x is normalized by the bath 
temperature Ti. The lines correspond (from above to below) 
to t; = 0.2; 0.3; 0.4; 0.5; 0.6; 0.7; 0.8. 



the intensity of the driving force. (For these systems, of 
course, the time needed to reach the stationary nonequi- 
librium state diverges as u — > 0.) In a sense, v plays a role 
similar to t^ in aging systems: the smaller v, the "older" 
the system. In a driven system TTI is satisfied, but FDT is 
not, even in the limit of vanishing driving force {v —* 0) M. 
Let us consider the slope x'(C) — dx/dC of the curve 
x(C). According to Q, for small enough driving forces, 
the effective temperature T°^ can be expressed in terms 
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of this slope: 

^"<^'-'^-(^)"'- <'' 

In all known cases one has T°^{C) > T. 

We can define a time scale T{q,v) by means of the 
relation: 

CiTiq,v),v)=q. (8) 

If q is larger than a threshold value (/ea (called the Ed- 
wards-Anderson order parameter) one has 

lim T{q, v) < oo. 
11— >o 



This is equivalent to the following definition of qEA'- 



qEA '■= lim lim C(t,w) 



(9) 



On the other hand, if g < qea, the time t{v) diverges as 
V goes to zero. 

Let us thus consider a thermometric measurement, 
performed on a characteristic time scale r. We wish to 
compare it with T°^ (r) , where 



T''*'{t) 



1 



X'{q) 



(10) 



q=C{T,v) 



where we have introduced the cross correlation function 

(7i2(t):= lim (Oi(t + r)02(i))c ■ (12) 

Using linear response, one has 

Oiit) ^ Oibit) + a f dtRi{t~t')02bit'), (13) 



02{t) ^ 02h{t) + a dtR2{t^t')0i^[t'), (14) 

where Ofa are the observables in the absence of coupling 
and Ri the corresponding response function (i — 1,2). 
At stationarity, to first order in the coupling a, the cross 
correlation function (n2h is given by 

Cl2{r)=af dT'Rl{T-T')C2{T') 



+ a dT'R2{~T')Ci{T' -t), (15) 



where we have introduced the correlation functions of the 
bare systems: 



Q(r):= lim(0,b(t + T)0,b(i))c 

t — >oo 

Then the rate of heat transfer writes: 



(16) 



3 Measurement procedure 

The measurement procedure is similar to the one described 
in Appendix C of B. It does not crucially depend on 
the nature of the thermometer, as long as it satisfies the 
fluctuation-dissipation theorem and has a tunable response 
time. We use a small but macroscopic thermometer in con- 
tact with a thermal bath at temperature T2. The driven 
system whose effective temperature is to be measured is 
in contact with a bath at temperature Ti. 

The thermometer is coupled to the observable Oi(5i) 
of the system via its observable 02{S2)- The interaction 
Hamiltonian writes: Hint ~ — a Oi (51)02(5*2). We remark 
that Oi and O2 are conjugate to each other. After a certain 
time, which depends on the coupling constant a and on 
l^*! — T2I, a stationary regime appears. 

We define as r™°'*'' the value of the temperature T2 
for which the net energy transfer between the system and 
the thermometer vanishes at stationarity. This tempera- 
ture is compared with the effective temperature T°^ of the 
system, defined by eq. (0). 

During the measurement procedure both the system 
and the thermometer are kept in contact with their own 
baths. The resulting system has two temperatures. (A sim- 
ple system with two temperatures has been introduced 
and discussed in ref. 0.) 

At stationarity, the net power gain for the thermome- 
ter, Q2, writes: 



Q2 = a {O1O2) = lim adrCi2{T), 



(11) 



dr (i?2(r)a,Ci(r) - i?i(T)a,C2(r)) . (17) 



From eq. (ph, we can substitute q for t as the integration 
variable: dq = Ci{T)dT. One can now exploit the fluctu- 
ation dissipation relations for the bare systems, namely 
Ri{t) = -C'i(T)/r<=ff(T) and i?2(T) = -C'2(r)/r2, and 
obtain 



dqR2[q) 



T2 



T-^iq) 



(18) 



The measured temperature T" 
which makes Q2 to vanish: 



is defined as the one 



T' 



s(^^)-l .^ JodqR2{q,q2)T^''{q)-' 
/p dTR2{q,q2) 



(19) 



We have introduced the parameter 52 := C'i(r2), where 
T2 represents the tunable characteristic time of the ther- 
mometer. T"^'^^{q2)~^ is the average of T'^'^(g)~^ weighted 
by R2{q, 92)- We remark that the measured temperature is 
independent of the coupling constant a, provided that it is 
small enough to ensure the validity of the linear response 
theory. On the other hand, T™°'*'' depends strongly on q2, 
because the lower boundary of the integrals appearing in 
eq. (19) is effectively cutoff at (72- 

For a system at equilibrium T'^^{q) = Ti for any q. 
Thus T'"°'''' — Ti as expected, whatever the characteristic 
time of the thermometer. 



Raphael Exartier, Luca Peliti: Measuring effective temperatures 



Let us consider a simple system with only two time 
sectors: 



ycff(^) 



Tu 



for T < Tea; 



Ti'>Ti, for T> TEA. 



(20) 



We have introduced the notation tea defined by the 
relation: Ci(tea) = 9ea- For g2 > Q'eAj which corresponds 
to probing the short time behaviour of the aging system, 
the lower boundary cutoff at g2 of the integrals of eq. ( p^ 
implies that the effective temperature of the driven system 
is constantly T'^^{q) = Ti over the integration interval and 
thus: 

T"^<=^''(g2)=Ti. (21) 

For 92 ^ (Zea the temperature T"''°^^ is not equal to T°^ — 
T[. Splitting the numerator of (nS) in two integrals from 
to gEA and from gEA to 1 we obtain 



rpW 



\Q: 



2) = 



J_Jl^dqlMg^ 
Ti J^ dTR2{q,q2) 
I 1 il^^dqR2{q,q2) 
Tl J^dTR2iq,q2) 



(22) 



The measured inverse temperature is a weighted average 
of the inverse temperature of the bath 1/Ti and the inverse 
effective temperature 1/T°^. Its value is intermediate be- 
tween them, as shown in figure 0, where the aging system 
and the thermometer are the ones described in the next 
section. 



an increasing function of q, we expect the measured tem- 
perature to be lower than the effective temperature when 
the thermometer probes the long time scales correspond- 
ing to the limit q -^ 0. Nevertheless, for intermediate time 
scales, if i?2(g, 52) is not peaked sharply enough around 
52, it is possible to measure a temperature higher than 
the effective one, as shown, e.g., in figure O . 
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Fig. 3. Comparison of the effective and measured temper- 
atures of the asymmetric spherical SK model for Ti — 10, 
J = 20 and v = .9. The solid line corresponds to the effective 
temperature of the bare model. The dashed line corresponds 
to the termometer measure. 
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Fig. 2. Comparison of the effective and measured temper- 
atures of the asymmetric spherical SK model for Ti = 10, 
J = 20 and v — .1. The solid line corresponds to the effective 
temperature of the bare model. The dashed line corresponds 
to the termometer measure. 



If we consider a driven system with many times scales, 
the measured temperature over a certain time scale 52 is 
the weighted average of all the effective temperatures of 
the system over this time scale. Since, on one hand, R2 is 
a decreasing function of q and, on the other hand, T is 



4 Effective temperature of a spherical SK 
model with randomly asymmetric bonds 

The previous considerations can be made more explicit in 
an exactly solvable model of a system-thermometer com- 
plex. 

We consider a spherical SK model with randomly asym- 
metric bonds [[10|. The Hamiltonian has the form 



1 ^ 
Hi := -- 2_^ JijSlSl + — 



2J = 1 



N 

E 



isi) 



(23) 



The "spin" variables SI can take any real value. The pa- 
rameter ri is a Lagrange multiplier which enforces the 
spherical constraint J2i=ii'^i)'^ — ^- Fhe interaction ma- 
trix Jf ■ is a symmetric matrix whose diagonal elements 
vanish and whose off-diagonal elements, for each pair {i, j} 
of indices, are independent Gaussian variables with zero 
mean and the following variance: 



ij^y 



.p 



1 



N 1 



(24) 



The parameter u > is a measure of the strength of the 
driving force (see below) . This spherical SK model is cou- 
pled to a single-spin paramagnetic thermometer via a bi- 
linear interaction of strength a. The Hamiltonian of the 
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paraniagnet is given by 
Ho 



r2 



{S2? 



(25) 



{dt+n{t))Ri2{t,t') =ai?22(i,t') 



The response time scale of the paramagnet is given by 
T2 := l/r2. Indeed, we recall that the bare response func- 
tion of a paramagnet has the expression 



i?2(i) = 6'(i)exp-''^*. 
The total Hamiltonian writes 

H := iJi(5i) + H2iS2) + i/i„t(^i, ^2), 
where the interaction Hamiltonian Hi^t is given by 

N 



Hi, 



:= —a 52 2_^ ^1- 



(26) 
(27) 

(28) 



Stability requires a < rir2. The dynamics of the system 
is described by a system of linear Langevin equations: 



dt"Rii{t,t")Rii{t",t'), 
(36) 
dt"Rix{t,t")Ri2{t",t'), 



(dt + r2)i?22(i, t') - aRuit, t') + 5{t - t'), (37) 

idt+r2)R2i{t,t') -ai?ii(i,t'). (38) 

The time ig can be freely chosen between the time the 
interaction was switched on and the observation times t 
and t'. The equations for the correlation function involve 
the response: 

{dt+n)Cn{t,t')=aC2i{t,t') + 2nRn{t',t) (39) 

+ [ dt"{j'^Rn{t,t")Cn{t",t') 

+ J^Cii{t,t")Rn{t',t")), 
idt+n)Ci2it,t')=aC22it,t') + 2TiR2iit',t) (40) 



dtSm^-—+b,{S,] 



vlit), 



dtS2{t) 



dH 

"dS'2 



■m{t)- 



In eq. (|29|), the driving field bi{Si) is given by: 



(29) 
(30) 

(31) 



+ 



dt"{J'^Riiit,t")Ci2it",t') 
+ J^Cn{t,t")R2i{t',t")), 



(dt + r2)C22{t, t') = aCi2{t, t') + 2T2R22{t' , t), (41) 
{dt + r2)C2i{t,t') ^ aCii{t,t') + 2T2Ri2{t\t). (42) 

After some time, the system enters a stationary regime 
where Qj{t,t') = C^j{t - t') and R-i.j{t,t') = R,j{t - t'). 
Choosing t^, t' and t in this regime, and taking advantage 
where J^^ is an antisymmetric matrix whose off-diagonal of Fourier analysis, one can solve the system for the re- 
elements, for each pair {i,j} of indices, are independent sponse and then the one for the correlation. For simplicity 
Gaussian random variables of zero mean and variance we suppose that the interaction between the system and 
equal to that of Jf, (eq. (gj)). The r]\ are thermal noises the thermometer has been switched on an infinite time 
at temperature Ti with zero mean and variance given by in the past. This corresponds to sending to — > —00. All 
{r]i{t)i]l{t')) — 2TiS{t — t'), while 772 is a thermal noise at the quantities depend on the value of ri which is cho- 
temperature T2. sen in order to verify the spherical condition: Cii{t,t) = 

In the thermodynamical limit {N -^ 00), it is possible j. ^^^ sf = 1. Taking the derivative of this condition 
to average over the disorder by the means of dynamical ^ith respect to time we obtain 
functional integration techniques. Thus the equations for 

the asymmetric spherical SK model in ( [29| ) reduce to a i;_ dCii{t,t) ^ ^.^ dCii{t,t) 

single equation for a single spin Si. The new system of 
equations reads 



lim 



lim 



0. 



(43) 



dtSiit) ^ -n{t)Si{t) + aS2{t) 



dt t'~,t+ dt 

Substituting eq. ( p9[ ) we obtain an equation for ri: 

n=T + aC2i{0) (44) 

+ (J2 + J'2) f dt" Rn{t,t")Cn{t'\t'). 

Jto 

Since C21, Rii and Cn all depend on ri this is an equation 
where J' = ^^(1 - v^)/{l + v^) J, and 771 is a renormal- for ri. We solved it numerically and then substituted the 
ized Gaussian thermal noise with zero mean and variance value of ri into the correlation and response functions. 



+ J'- 



dt'Rn{t,t')Si{t')+m{t), 



dtS2{t) = -r2S2{t)+aSi{t)+f]2it), 



(32) 



(33) 



After this step the solution is completed and we can search 
for the temperature T2 of the paramagnetic thermometer 
(34) which makes the heat flux to vanish. 

The power exchanged between the thermometer and 
For the response we obtain the following autonomous the system at stationarity is given by 
system: 



given by 

im {t)m {t')) = 2TiS{t - 1') + j^Cii (t, t'). 



{dt + ri {t))Rn{t, t') = aR2i{t, t') + 5{t - t') 



(35) 



Q2- = 



dHn 

dSi 



Si)- 



dS2 



S, 



(45) 
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From the definition (E8h of ifint we obtain 

Q2(t) = -ab2(t) — — h aSi{t) -— 

at at 

= lim {~adtC2i{t,t')+adtCi2{t,t')) 



(46) 



By using the resuhs of the appendix and taking an inverse 
Fourier transformation, we finally obtain 



— iuj [Cuiuj) - C2i(a;) 
dio 



-2a 



2tt 



■ Lo Im 



(6-12 (^)) 



(47) 



The temperature measured by the thermometer is the 
one which makes the heat flux Q2 to vanish. The re- 
sulting measured temperatures are shown in figure ^ and 
should be compared with the expected effective tempera- 
ture shown in figure H. 



5 Results 

Figures ^ and |5| show the effective and measured values of 
the temperature as a function of the value of the correla- 
tion function, for the model under study. The two quan- 
titites behave similarly, in that they are close to the equi- 
librium value Ti for q > qea, and start increasing, as q 
becomes smaller and smaller, for q < qea- Nevertheless 
it is possible to identify a quantitative discrepancy, anal- 
ogous to the one discussed in section 2. For small values 
of the asymmetry parameter v, the driven system exhibits 
two clearly separated regimes, with T°^ — Ti for q > gEA, 
and a temperature increase for q < qEA- For high values 
of q, T"^'^'^ remains close to Ti, but is much smaller than 
T°^ in the low-g region, as shown in fig. 0. For higher 
V, T°^ increases smoothly all along the range of q. Being 
an average of T'^^ over a given time scale, T"^'^'^^ becomes 
quickly higher than Ti as it feels the increasing of r°*^ 
even for short time scales. Then it keeps on increasing, 
at a slower pace than T"^ , as shown in fig. pi Figures ra 
show the ratio ymcas^ycff g^ ^ function of the asymmetry 
parameter v. The ratio remains close to the ideal value 
1 only for q > q^A- Ideally, similar plots would apply to 
an aging system as a function of the inverse waiting time, 
because of the correspondence between v and the waiting 
time discussed above. 

In figure |7| and |8| the behavior of the measured and 
of the effective temperature is shown as a function of the 
asymmetry parameter v is shown for different values of the 
thermometer response time T2. The figures are in qualita- 
tive agreement only for short time scales which are rep- 
resented by the lower lines of the plots. When the ther- 
mometer probes longer times scales, T"^°'^^{t2) shows a 
non monotonic behaviour which does not appear in the 
T'^^{t) plot. In the waiting time representation, a ther- 
mometer with a fixed (but long) reaction time T2 would 
first yield higher and higher readings, as it attempts to 
approach the flat part of the graph shown in fig. |l|, but 
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Fig. 4. Effective temperature T'^ vs. q for the bare asym- 
metric spherical SK model for Ti = 10, J = 20, and for dif- 
ferent values of the asymmetry parameter v. T^ is normal- 
ized by the temperature Ti of the thermal bath coupled to 
the aging system. The lines correspond (from above to below) 
to w = 0.1; 0.15; 0.2; 0.3; 0.4; 0.5; 0.6; 0.7; 0.8; 0.9. The Edwards- 
Anderson order parameter qsA of the corresponding symmetric 
model is given by qea — Ti/J — .5, and corresponds to the 
value of q from where the curves diverge. 
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Fig. 5. Measured effective temperature r™"=^= ys. q of 
the asymmetric spherical SK model for Ti = 10, J = 
20, a — .1, and for different values of the asymmetry 
V. The lines correspond (from above to below) to t; = 
0.1; 0.15; 0.2; 0.3; 0.4; 0.5; 0.6; 0.7; 0.8; 0.9. 



will eventually read the temperature of the thermal bath 
as tyj :^ T2. This leads to a non monotonic behavior as a 
function of t^, and an analogous one as a function of v, 
at least for sufficiently slow thermometers. 

In conclusion we have shown in an exactly solvable 
model of a thermometric measurement in a "glassy" sys- 
tem that, while thermometer do indeed measure temper- 
atures higher than the one of the environment if they are 
slow enough, the relation of the measured temperature 
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with the effective temperature defined, e.g., in [H is far 
from trivial. 
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Fig. 6. Ratio r'^'=^Vr° 



asymmetry parameter v of the 



asymmetric spherical SK model for T\ = 10, J = 20, and for 
different values of the overlap q. The lines correspond (from 
above to below) to q = 0.8; 0.5; 0.4; 0.1; 0.01; 0.001. 



jmeas/J 




Fig. 7. Measured effective temperature T™'='"^ vs. asymmetry 
parameter v of the asymmetric spherical SK model for T\ = 10, 
J = 20 , a = .1, and for different values of the parameter T2 
which sets the characteristic time of the thermometer. r™<=^= 
is normalized by the temperature T\ of the thermal bath cou- 
pled to the aging system. The lines correspond (from above to 
below) to ra = 100; 10; 1.0; 0.1; 0.01. 
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Fig. 8. Effective temperature T'^ vs. asymmetry parameter v 
of the asymmetric spherical SK model for Ti = 10, J = 20 and 
for different values of r. T° is normalized by the temperature 
Ti of the thermal bath coupled to the aging system. The lines 
correspond (from above to below) to r = 100; 10.; 1.; 0.1; 0.01. 
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Appendix 

Taking the Fourier transform of the autonomous system 
for response functions, equations (p5[-p8|), we obtain 

{iuj + ri)Rii{uj) = aR2i{uj) + f^Rniuf + 1, (48) 
{luj + ri)Ri2{oj) = aR22{uj) + J'^Rii{oj)Ri2{ijj), (49) 
{iuj + r2)R22{^)^aRi2{uj) + l, (50) 

{iuj + r2)i?2i('^) = aRii{uj). (51) 



Again, from (p9H42) we obtain the following equation for 
the correlation functions: 

{ioj + ri)C'ii(cj) = aC2i{i^) + J'^Rii{uj)Cii{Lo) (52) 

+J2(7ii(w)i?ii(cj) + 2TiRn{uj), 
{iuj + ri)Ci2{uj) = aC22{i^) + J''^Rii{uj)C\2{uj) (53) 
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+J^Cii(a;)i?2i(w) + 2T^R2i{u;), 
{luj + r2)C22{uj) = aC'i2(w) + 2r2i?22(w), (54) 

{iuo + r2)C2i{t^) = aCii[uj) + 2T^R^). (55) 

In Fourier space, the equation (^J) for the spherical pa- 
rameter Ti writes: 



/T 
2^C2i(w) 



+ (J2 + J'2) 



dw - 



27r 



i?ii(a;)Cii(tj). (56) 



Equations (E8-56) form a nine-equation system for the 
response and correlation functions and for the spherical 
parameter, which is possible to solve explicitly. Equations 
(p8| ) and ( ^ ) yield a second-order algebraic equation for 

i?ii(w): 



J'^R{^ ~ [iuj + ri 



ito + r2 



Rii + l = 0. 



(57) 



We choose the solution which respects the symmetries of 
Im(^ii(u;)), Re(^ii(a;)), and recovers the right value for 
i?ii(a;) in the limit a ^ 0. 

From equations (Eoh and (Eo) we obtain 



^12(0;) = -R2l(w) 

while equation (^^ yields: 
1 



iuj + r2 



i?ii(w). 



i?22(w) 



iuj + r2 {ioj + r2)^ 
From equations (p2|) and (B3) we obtain 



Cii = 

C21 = 



2Ti + 2T2aV(a;^ + r|) 
|i?ii|-^-J2 



(r2 - iLi;)Cii + 2r2^ii(^) 



(58) 



Riiito). (59) 



(60) 
(61) 



As shown by equation (|47|) , the measurement procedure 
of T"^'''*^ gives a particular relevance to the imaginary part 

of (721 (w): 

T (r^ \ 2r2lm(i?n)+c^Cii .„„. 

Im(C2i) = -a 2- — 2 (62) 



